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ABSTRACT 

A statistical theory of inelastic neutron energy spectra for intermediate 
excitation energies is developed. The final state of the residual nucleus may 
be either discrete or in the continuum. Calculated neutron energy distributions 
for aluminum are presented. 
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DEFINITION OF SYMBOLS 


Primes 

c 


a 


P 

E 

P 

E * 
P 

S 

P 

Q 

I 


i 

s 

F 

j 
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refer to the final state. 

The set of quantum numbers {a, I, i, s, I , J, n } which 
specifies the reaction channel. 

The set {p, E^} which labels a pair of particles and their 
state of excitation. 

The pair of particles. 

The eenter-of-mass kinetic energy of the pair p. 

The energy of excitation of the pair p. 

The separation energy of the pair p from the compound system. 
The energy released by the reaction. 

The intrinsic spin of the target. All spins are in units of Ti , 
Planck's constant divided by 2 7 r . 

The intrinsic spin of the bombarding particle. 

Toe channel spin; s = I + i . 

The orbital angular momentum of the pair. 

The total angular momentum; J = l + s = jf + I + i. 

The z component of J. 

“ft 

The total parity (the product of ( 1) with the intrinsic parities 
of the pair a ) . 


G , Let the decay and production of compound systems for fixed 

J and r be equal, so that at any time there are N systems 

present. Then NG^/fi is the number of decays per second 

cx 


k 


The wave number for the relative motion of the pair a . 


P c < E*,I) 


P d (E*,I) 


The number of suites in the continuum with angular momentum 
I and excitation energy in dE* at E* . 

The level density for discrete suites. 


__ c 

a 

fV fY 


The average cross section with the residual nucleus excited 
into the continuum region. 


_ d 

a 

a rv 


H(X) 

K 

C’ 

g n 


The average cross section for a discrete final suite. 

If E < E the residual nucleus will be in a discrete state 
P c 

If E > E the residual nucleus will be excited into the 
P c 

continuum region. 

For any X if X < 0 then H(X) 0. If X > 0 then H(X) 1 
The moment of inertia of the nucleus. 

A constant depending upon the given nucleus. 

The density of neutron orbits at the Fermi level. 

The density of proton orbits at the Fermi level. 


G 
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The total density of orbits at the Fermi level. 

The nuclear ’temperature,” not to be confused with the 
’’temperature" in the simple evaporation model. 


The Fermi energy; the excitation energy above the completely 
degenerate state. 


J. 


i Z 


The atomic mass number. 


.th 


The Z - component of angular momentum of the i particle 
or hole. 






<X> 

The 

f(U ) 

V 

The 

*<V 

The 

m 

The 

V 

The 

W 

The 

i 

The 

V 2 

The 

D(U) 

The 

gas 


average value of X. 
probability density function of U 
probability density function of J . 

Tj 

reduced mass of the pair 

real part of the optical model potential 

imaginary part of the optical model |X)tential 

total wave function 

Laplacian 

nuclear energy level density according to the Fermi 
model 


INTRODUCTION 


Space Vehicle Radiation Shielding 

Space vehicles are exposed to proton fluxes produced by solar flares, 
galactic cosmic sources, and the trapped radiation belts 1 1 J . Proton collisions 
with the vehicle produce secondary neutrons which are then involved in sub- 
sequent nuclear reactions. At neutron energies between 8 and 18 MeV about 
one half of these secondary neutrons will initiate inelastic neutron scattering 
reactions . The energy spectra of these inelastically scattered neutrons is 
required ls part of the input for Monte Carlo calculations which predict the 
radiation dose rate inside the space vehicle. 

Nuclear Reaction Mechanisms for Inelastic Scattering 

The first problem of nuclear reaction analysis is to establish the 
mechanism of the interaction. The reaction mechanism is determined by the 
number of particles which share the energy introduced into the target nucleus. 
Nuclear reactions are classified for convenience as direct, intermediate, and 
statistical reactions. 

If only a few degrees of freedom of the target nucleus are excited, then 
the reaction is direct ( 2 ] . The momentum transferred to the residual nucleus 
is small. This commonly results in angular distributions with large asym- 
metric peaks in either the forward or the back angles. 

Doorway or intermediate mechanism states are iormed as additional 
nucleons share the energy [3j. The doorway states may either collapse 
rapidly to single-particle states or proceed to more complicated configurations. 
Reactions which go to completion via energy excitations of intermediate com- 
plexity have not received adequate theoretical treatment. 

If the complexity of the intermediate state continues to increase, then 
ultimately all nucleons in the nucleus share the available energy. Thus a 
compound nucleus is formed |4] . The angular distributions resulting from 
compound nucleus decay are symmetric about 90° in the center-of-mass system. 

At present only two of these three types are available for treating 
inelastic scattering: the direct reaction theory and the compound nucleus or 
statistical theory. In this report it is assumed that the inelastic neutron 
scattering proceeds via compound nucleus formation. 


THE COMPOUND NUCLEUS 


A compound nucleus is a metastable system which is formed when an 
incident particle is captured by a target nucleus, and the available energy is 
shared by all nucleons in the system. If one makes the Bohr |5J assumption 
that the compound nucleus is analogous with an equilibrium system in classical 
thermodynamics , then the decay of the compound nucleus should be dependent 
of the mode of creation. This assumption leads to a factorization of the 
reaction cross section into two independent elements (assuming a fixjd total 
angular momentum) . The first of these is the cross section for the formation 
of the metastable compound system, anti the second is the probability for sub- 
sequent decay through the channel of interest. In this report only the exit 
channel containing the inelastically scattered neutron is considered. 

If no external magnetic fields are present, a time-reversed state must 
also satisfy the equations of motion which govern the decay of the compound 
system. This symmetry is expressed mathematically by the reciprocity 
theorem [4j. The rcjipiocity theorem in conjunction with the Bohr assump- 
tion allows the compound nucleus decay probability to be expressed in terms 
of the compound nucleus formation cross section for the inverse reaction. 

In the inverse reaction the target is in an excited state arid there is no way to 
experimentally measure the inverse compound nucleus formation cro3S section. 
It is assumed that this cross section is the same as the compound nucleus for- 
mation cross section for a target nucleus in the ground suite. 

These assumptions allow’ the inelastic neutron cross section to be 
expressed in terms of the initial state compound nucleus formation cross 
section. The compound nucleus formation cross section is approximated by 
the optical model [6] reaction cross section. 

For a given excitation energy several states in the compound nucleus 
will contribute to the total reaction amplitude. The statistical hypothesis [ 7 { 
assumes that the interference terms between the different transition amplitudes 
cancel. 


The compound nucleus decay probability also contains a term which 
represents the available final -state phase space. It is normally assumed that 
the density of final states of the residual nucleus depends upon the excitation 
energy and spin, but is independent of the parity of the state. The energy and 
spin dependence are usually separated into independent multiplicative factors. 


o 



The energy dependence is approximated b> using the Fermi gas model of the 
nucleus (8), and the smn dependence for multiparticie excited suites is approxi- 
mated by a Gaussian distribution |!)j. 

Thus the inelastic neutron cross section is expressed in terms of func- 
tions which are determined by optical model and level density parameters. 

The physical implications and validity of the above assumptions will be dis- 
cussed in the sections which follow. 


THE HAUSER-FESHBACH MODEL 


Early evaporation models which were developed by liethe and Weisskopf 
considered only conservation of energy. The Hauser-Feshbach 1 10] model 
presented below includes conservation of total angular momentum and parity. 

The incident neutron energies are spread over some small energy 
interval. This implies excitation of compound states within a corresponding 
energy interval. The statistical assumption states that there are many excited 
states in this interval, a»iu the widths and energies of the excited states are 
randomly distributed within this energy interval. The cross sect* n calculated 
is the average over individual resonances in this energy interval. 

Let <7^^, be the average angle -integrated cross section for transition 

from state a to state a' . This cross section is the sum over all partial 

Jtt 

cross sections, ^ , , which conserve angular momentum, J, and parity, 7 r . 
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It is assumed that each partial cross section may be factored into two 
factors. The first factor is the cross section for the formation of the com- 
pound nucleus, ctJo) , and depends only upon the initial state. The second 


factor is the branching ratio, G V G^ f , and depends only upon the final 
state. 01 ' Q 


( 2 ) 
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cr , ill a G /' G , f 
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The reciprocity theorem may be used to express the branching ratio 
in terms of the compound -nucleus -formation cross section. 


2 -Jr 

C7 

a a a 


k 2 , a J7 [ 
a q a 


( 3 ) 


From equations (2) and (3) one obtains 


J,r v ~Jr /V ~ j7r i 2 J,r v ~«*r / \ ' ^ j7r 
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Equating terms in a gives 


.Jr 

i 
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k 2 a ^(a) 
a c 


(4) 


Each a (<*) is obtained from the optica 1-model absorption cross 
c 


section, <j 


a a 


_ \ J ,ir . . 

a n * A (fY ) ( 5 ) 

a a t - c 
Jr 

In terms of the optical -model complex phase shifts, 6^ ^ this is 

a * (r/k 2 ) \ (21 + 1) ( 1 - I exp 2 i 6 I 2 ) , (G) 

a a a at 

t 

where f is the orbital angular momentum of the interacting particles. For 
the partial cross sections one has 


J,r 

o c (a) 


( r (2J + 1) 


k* (21 + 1) (2i 1) I 
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14 i 

X 

S = I i - i| i = jj-s j 


J 4 S 

T (1 - Iexp2i 6 I 2 ) 


at 
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The transmission functions are defined by 


T. <rv) * 1 - I exp 2 i 6 A 2 . 
i at l 


<«) 


Thus, from equations (4), (7), and ( b) , one finds 

C j j |»(2J+ 1)/(2I+ 1)(21 + 1)J V (a) 

<,s 

From equation ( 2) , 

a Jlr , = |tt( 2J + l)/k 2 (21 + 1) (21 + 1)| V T (a) 
at a a . Lj i 

From equations (1), (9), and (10), 
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Since it is not possible to define the final energy more accurately than 
the initial energy, there is an energy interval associated with the final state. 

If the compound system is highly excited, then the residual nucleus is likely to 
be highly excited also. For a highly excited residual nucleus the level density 
is represented by a continuous density function. However, if the compound 
nucleus has a lower excitation energy, then the residual nucleus may be in a 
discrete energy state. In general there will be some transitions to continuous 
final states and some transitions to discrete final states. The experimentally 
measured neutron spectrum is the average cross section times the number of 
final states in the energy interval. Thus the energy spectrum of neutrons 
emitted from the compound nucleus is 


N(E .) = pJE* ,I')5 C , H(E - E ,) + p,(E* ,I-)5 d ,H(E -E), 

p ^ a at' c p’ d eta' p’ c 


( 12 ) 


where H is the Heaviside function. 
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The discrete density of states, p^ , is unity if dE is chosen so small 
that only one final discrete state is excited. The continuous density of states, 


p , is taken to be the form 

c 


( E* , I) = p(E*)(2l + 1) exp|- 1(1 + 1)/2 ij 2 | 


( 13) 


The above spin level -density formula predicts too many levels of high 
spin at low energies. All levels of high spin at low energies should be excluded 
if they are not energetically obtainable by adding rotation energy to the ground 
state. If this spin cutoff is included the neutron energy spectrum becomes 
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where g = ( 2J + 1) /( 21 + 1) ( 2i + 1) 
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NUCLEAR LEVEL DENSITY 


The inelastic neutron energy spectrum contains a factor which represents 
the available final -state phase space. For a two-body reaction, the phase space 
factor is determined predominately by the level density in the residual nucleus. 

Nuclear level densities exhibit an essentially exponential increase with 
excitation energy. This rapid increase may be basically explained as a result 
of the addition of many elementary one-particle excitations to form a highly 
excited many-particle level. This is a characteristic of the Fermi gas model 
which only required that the excitation energy of the level be distributed among 
several excited nucleons in a manner consistant with the Pauli principle. 

Nuclear level densities are strongly influenced by nuclear shell struc- 
ture. In particular, for a fixed excitation energy, the level densities of magic 
nuclei are orders of magnitude smaller than the level densities of adjacent 
nuclei. 


Odd-even effects which are attributed to differences in pairing energies 
are also present. At low excitation energies, the odd nuclei often exhibit 
several shell model levels, but even nuclei have few levels in this region of 
excitation. Above the region of collective excitations, the even nuclei have 
many levels which arise from the numerous ways in vvhich unpaired nucleons 
can be arranged. 

The above effects are most pronounced at low energies and decrease 
with increasing excitation. However, these factors do not entirely disappear 
and must be considered in any satisfactory theory of nuclear level densities. 


Energy Distribution 

Models which describe nuclear level densities must be related to Loth 
the Fermi gas model and to the shell model. The odd -even effect indicates an 
additional complexity. Newton [8] has considered the nucleus as a Fermi gas 
and included shell effects by replacing the density of nucleon orbits at the 
Fermi level by the density of nucleon orbits in the shell model. The equivalent 
energy of the Fermi gas is obtained by subtracting the pairing energy correc- 
tion from the nuclear excitation energy. In this model the level density is 
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D ( U) = C'AG 2 G n 2 (2U ♦ 30 ) 2 exp | -2( ir 2 G U/6) /2 J , 


where 0 


m 


/2 


and G = G VT + G 
N z 


( 15) 


Effective values of G for various nuclei have been tabulated by Cameron [ 1 1 j . 


Angular Momentum Distribution 


Consider a Fermi gas with v excited holes and particles. The projec- 
tion of the total angular momentum on the Z-axis of the system is 
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i' 

V 

_ * 
i=l 


J. 


1Z 


(16) 


The random variables J , . . . , J are independent and have the same dis 

lz vz 

tribution with standard deviation o i- 0 . 


< J iz > ■ 0 
<hz> - 


(17) 


Define the arithmetic mean of the v independent random variables as 


U = J /v . 
V z 


(1H) 


It immediately follows from the central limit theorem [12J that U has an 
asymptotically normal distribution with mean ^ and standard deviation 


(J fSV . 


The normal distribution is defined by 


f(x) = 


(x - <x> ) 


2 -i 


\/ 27T <x 2 > 


exp 


2 <x 2 > J 


(19) 
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and 
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f(U ) = 

si 2ir a 2 /v 


exp 


u- 


2 a / v 
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( 20 v 


Since f( U ) dlJ = g(J_)dJ it follows that 

V V A A 
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f(U ) 
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si 2tt V (T 4 


exp 
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2vg 

L v J 


( 21 ) 


The density of states of spin I, p(I) , is the difference between the level 
density of spin J = land J = I + 1 . 

A La 


9 (I) 


g(I) - g(I + 1) & 


(21 + 1) exp[- 1(1+ 1)/2 ct 2 ] 
2(2tt) /z a 3 


( 22 ) 


The parameter a 2 



is called a spin cutoff parameter, which must be 


determined for a given nucleus. The spin cutoff is a function of the Fermi 
energy. However, a is normally taken to be constant, since if it is energy 
dependent, the energy integrals and spin sums in the neutron inelastic energy 
spectrum cannot be separated. 


If p(I) is energy independent, then it predicts the same ratio of high- 
spin states and low -spin states at ail excitation energies. This is incorrect 
since high spin states do not appear at low' energies. The low-lying levels re 
shell model states of low r angular momentum. Rotational bands are built upon 
these low -lying states, and the rotational energy determines the high-soin 
states. The rotational energy is given by 

E = ( Ti 2 /2K) 1(1+1) . (23) 


If this rotational energy is added to the ground state, then determines the 

energy level at which a high -spin state first appears. This fact is represented 
in the neutron inelastic energy spectrum by an energy cutoff which prohibits high 
spins at low energies. 
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Parity Distribution 


Single Fermion states have parity ( + 1) or (-1) . The resulting parity 
of the system may be determined from elementary probability considerations. 

It may be shown that only a small number of negative parity states or a small 
number of positive parity states are required to obtain a nearly equal probability 
for positive and negative total parity. 

THE OPTICAL MODEL 


The optical model represents the interaction of a nucleon with the nucleus 
in terms of a single-particle complex potential. The real part of the potential 
describes the scattering of nucleons; the imaginary part describes absorption. 

The Schroedinger equation for neutrons can be written as 


v 2 ii + - (V + iW) ] ip = 0. 

The solution is obtained by expanding if, in terms of Legendre polynomials. 

u = X C,(r) PI coso) . 

( 

Then 4 ( \r) satisfies the equation 


1 

7 dr 


dr 


+ -fr- [E - <v + iW) j ^ - - ( -- yr ■ 1 * = o . 


The general solution of this equation which asymptotically approaches the 
equation of a free particle is of the form 




= c i 
i a 


C2 ^ 
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The boundary condition at the origin is 

lim [ r ^ ( r) ] = 0 . 
r— 0 
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This condition determines the complex phase shifts 6^ f which are defined by 


£jL 

C 2 


= exp(2i 6 ) £ n = lim 

I — >0 


- r* 


II 




The asymptotic expansion is 


h = C2 


!(*,“ * *,') 


- ( 1 - exp ( 2i 6 ))^y . 

at * | 


t u 


The first term represents the S. * term in the incident plane wave expansion. 

th 

The second term is the f -order scattered wave. 


The reaction cross section, a , is the ratio N /N, where N is the net 

r a a 

flux into a large sphere of radius r 0 as computed from the complete wave 
function, and N is the incident plane wave flux. 


o = 

r 


2imV 


/ (“l^ ~ ro sino do d 0 . 


The cross section calculated by the optical model is an average cross 
section over an energy region containing many resonance levels. The optical 
model is satisfactory when the width of the beam is considerably larger than 
the level spacing. Thus, over an energy interval containing many resonances. 


/ , v. dE 
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The average reaction cross section is 
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r k 


7T 


V ( 2£ + 1) ( 1 - IT?) 

n l 


or 


a = , 

r k 


T " L (2f + 1) 


ot a 


a - h f i 2 > - ^ i * > 


li 


The second term results from resonance reactions for which differs 
significantly from n . The first term is the absorption cross section for com- 
pound nucleus formation, i.e., the absorption results in formation of a com- 
pound system with excitation energy in a region which contains many resonance 
levels. Thus 

a a a = 1)(1 - 'V 2) • 

Ot i 

Various functional forms have been used to approximate the optical 
model potential 6, 11. For incident energies between 5 and 10 MeV, it is 
best to work with diffuse boundary potentials. In this case, it is necessary to 
use a computer in order to 3olve the Schroedinger equation and obtain the phase 
shifts. 


CALCULATIONS 


The inelastic neutron cross section for aluminum for an incident energy 
of 5 MeV is compcsed primarily of transitions to discrete energy levels. The 
results of the computer program ABACUS [ 12 J for an incident energy of 5 MeV 
are presented in Table I. 

As the incident energy increases the transitions are to both discrete 
final states and final states in the continuum. Figures 1 and 2 show these 
mixed transitions at 9 MeV for 18 and 10 discrete levels, respectively. Figure 
3 is the result for an incident energy of 12 MeV. 

At higher incident energies the cross section is dominated by transi- 
tions to the continuum. Figures 4, 5, and 6 include continuum transitions only. 

Figures 1 through 6 are not the neutron energy spectra as represented 
by Equation ( 14) , but are a composite of ABACUS and a computer code 
developed by R. Snow and M. C. George ( 13 J . The curves are all normalized 
to an integrated cross section of 0. 8 barns. Some of the options available in 
the continuum code are indicated in Figure 7. 
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TABLE I. DISCRETE TRANSITIONS OF INELASTICA LLY SCATTERED 
NEUTRONS FOR AN INCIDENT ENERGY OF 5 MeV 


DISCRETE ENERGY 
LEVELS (MeV) 

CROSS SECTION 
(BARNS) 

0. 842 

0.053547 

1.013 

0. 109156 

2.212 

0. 173929 

2.731 

0. 116153 

2. 976 

0.064010 

3.010 

0. 187355 

3.674 

0.022831 

3. 951 

0.045302 

4.052 

0.017250 
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FIGURE 1. APPROXIMATELY CALCULATED CONTINUUM ENERGY DISTRIBUTION 
AND EIGHTEEN DISCRETE TRANSITIONS FOR INELASTIC A LLY 
SCATTERED NEUTRONS WITH AN INCIDENT ENERGY OF 9 MeV 
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0123456789 


E'(M*V) 

FIGURE 2. APPROXIMATELY CALCULATED CONTINUUM ENERGY DISTRIBUTION 
AND TEN DISCRETE TRANSITIONS FOR INE1.ASTICALLY SCATTERED 
£ NEUTRONS WITH AN INCIDENT ENERGY OF 9 MeV 
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LATELY CALCULATED CONTINUUM ENERGY DISTRIBUTION 
:RETE TRANSITIONS FOR PCELASTICALLY SCATTERED 
RONS W ITH AN INCIDENT ENERGY OF 12 MeV 
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FIGURE 4. EXACTLY CALCULATED ENERGY DISTRIBUTION OF INELASTIC ALLY 
SCATTERED NEUTRONS FOR AN INCIDENT ENERGY OF 15 MeV 






FIGURE 6. EXACTLY CALCULATED ENERGY DISTRIBUTION OF INELASTICALLY 
SCATTERED NEUTRONS FOR AN INCIDENT ENERGY OF 21 MeV 
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